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We extend our ensemble density functional approach to quantum Hall systems to include non- 
coUinear spins to study charge-spin textures in inhomogeneous quantum Hall systems. We have 
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factor as a function of the stiffness of an external confining potential. For soft enough edges, these 
systems reconstruct to a state in which the electron spins rotate gently as the edge is approached, 
with a non-trivial spin-charge texture at the edge of the system. 
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I. INTRODUCTION 



It has recently become cleatthat the spin degree of freedom plays a significant role near ferromagnetic fillings in the 
quantum Hall effect (QHE)ErB. This is because of the low ratio of the Zeeman energy Ez = g* ^irB to the Coulomb 
energy Ec = e'^/{€olB)- Here, g* is the Lande 5-factor, fiB the Bohr magneton, B the applied magnetic field strength, 
eo the static dielectric constant, and Ib ~ y^hc/eB the magnetic length. For GaAs systems, the low values of the 
Lande factor and of the electron effective mass m* conspire to make g — Ez/Ec « Ez/{fiLUc) ^ 0.02 for fields in the 
range of a few tesla to about ten tesla. Here, uJc = eB/{m*c) is the cyclotron frequency. Nonetheless, single-particle 
spin-flip excitations still cost a large amount of energy, because of the loss of exchange energy associated with a 
spin flip. This means that the spin degree of freedom is controlled by the inter-electron Coulomb energy, and not 
by the Zeeman energy. One consequence is that if a single spin is reversed, it becomes energetically favorable for 
the system to smoothly rotate the magnetization direction to restore it over some distance from the reversed spin. 
Due to the connection between flux and charge density in incompressible ferromagnetic QHE ground states, such 
spin textures also acquire a charge density, and the resulting spin-charge textures are commonly called 'skyrmions'. 
(More accurately, skyrmions are the particular type of spin-charge-textures that show up in cettjain continuum models, 
such as the nonlinear cr-modelcl.) There is now ample theoreticalnO'Q and experimental worluTa suggesting that such 
skyrmions are indeed the low-energy charged excitations, at least near filling factor v = 1— This is, for example, 
manifested in the rapid destruction of the ground state polarization observed experimentallyQ as the filling factor is 
varied away from unity. Recent theoretical workt9il3 has also indicated that the edge reconstruction of ferromagnetic 
QHE systems may acquire non-trivial spin textures (or charge density waves) as the edge confinement is softened. 
This may radically alter our understanding of edges in QHE systems, and, concomitantly, our interpretations of 
experiments which probe the gapless edge modes. This is at the present best understood for edges of v = 1 systems. 
One question we wish to address here is whether similar charge-spin textures can occur at the edges of systems 
in regimes of the ferromagnetic fractional quantum Hall effect (FQHE) (such as = 1/3), where not just electron 
exchange but also correlations are important. nn 

We have previously developed an ensemble density functional approach for spin-polarized systemsllj'tj, and subse- 
quently generalized that approach tojaclude the electron spin, but with the spin quantization axis constrained to be 
parallel to the external magnetic fieldll3. We present here a further extension which is able to deal with a rotating spin 
quantization axis. Advantages of our DFT approach are that it includes electron interactions beyond exchange, and 
that it can be applied to large inhomogeneous fractional QHE systems. This makes our ensemble spin DFT approach 
the only available method which can be applied to general inhomogeneous QHE systems, spanning regimes from the 
FQHE to the semiclassical, and which includes the spin degree of freedom and Landau-level mixing. We have used 
this approach to study the edge reconstruction of circularly symmetric quantum dots. Our results show that as the 
edge confinement is softened, the system goes from a spin-polarized sharp edge to a softer edge with a non-trivial spin 
texture, in agreement with_L-esults obtained by other groups using the Hartree-Fock approximation (at integer filling) 
or field-theoretical modelsll3. A new result here is to show how the reconstruction to a spin-charge textured edge can 
also happen for FQHE systems. 

This paper is organized as follows. In Section ||, we present the general ensemble density functional theory for 
non-coUinear spins. Section HI presents some technical details of the theory, including an extension of previously used 
exchange-correlation energies for QHE systems to include both higher Landau levels and electron spin. In Section 
[V we present results from numerical calculations of quantum dots at unit bulk filling factor and at 1/3 bulk filling 
factor. Finally, section ^ contains conclusions and a discussion. 



II. ENSEMBLE SPIN DENSITY FUNCTIONAL THEORY FOR NON-COLLINEAR SPINS 

In its most general formulation, spin density functional theory, as developed, by von Barth and Hedin0, allows for 
non-coUinear spins. This is based on a generalization of the Hohenberg-KohnE3 theorem in which the electron density 
n(r) is replaced by the single-particle density matrix ^^.^'(r) = {0\ijl{r)il:„i {r)\0) , where V'o-(r) (V'J.(r)) is the usual 
electron annihilation (creation) operator for an electron of spin a at position r, and |0) is the ground state. We will 
use notations in which a = ±1 or a =tii: with up-spin (t) corresponding to tr = —1. A computationally useful 
approach is then constructed in the usual way by considering an auxiliary non-interacting system in some effective 
external potential Vs{r) chosen so that this system has the same ground state single-particle density matrix p^ ,ji[r) 
as the interacting system at hand. A vpJational principle associated with the generalized Hohenberg-Kohn theorem 
then yields Kohn-Sham (KS) equationaij, which now include spin-dependent exchange-correlation potentials 

W(r)^^%%^. (1) 
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A difficulty is that one does not usually have reliable approximations for the exchange-correlation potentials l^co-o-' (r), 
not even in the local density approximation (LDA) . Consequently, the density matrix is usually assumed to be diagonal 
for all r, which means that the direction of the magnetization is assumed to be constant, and that direction is then 
conveniently chosen as the spin quantization axis. Thus, only up- and down-spin densities enter into the KS equations, 
and for the LDA (or extensions including generalized gradient approximations), one only needs to know the exchange- 
correlation energy Exdn,^) of a uniform system of density n ^ + and polarization ^ = (n^ + 7i|)/(rt| — 
There-Epdst now very accurate calculations of E^dn, ^) for two- and three-dimensional electron gases (in zero magnetic 
field)B0. 

However, the approximation of constant magnetization direction obviously does not work in systems fcp-which it_is 
known that the magnetization direction changes in space. Examples of such systems are MusSnEj, 7-Fecj, U3Pt4EJ, 
and QHE systems near unit filling. For such systems the full single-particle density matrix has to be used, and the 
problem then arises SAjta how one should construct a LDA. We will here follow an approach developed by Kiibler, 
Sticht and co-worker£3E3, and extend their approach to an ensemble DPT appropriate for QHE systems. The basic 
idea is to locally rotate the spin quantization axis to obtain a representation which locally diagonalizes the single- 
particle density matrix. The advantage of this procedure is that in order to construct a LDA, one then only needs 
the exchange-correlation energy as a function of spin-up and spin-down densities (or total density and polarization), 
for which approximations exist. The price one has to pay is to introduce local spin rotation angles 9{r) and 
which jComplicates the KS equations. This approach has given good results in applications to materials such as 7-FeE3, 
U3Pt40, and a-Fe203O. It can also give the spin stiffness important in studying spin-charge textures in the QHE. 
That this can happen in the LDA is not immediately obvious - in field-theoretical approaches the spin stiffness enters 



via a gradient of magnetization, and such gradients are neglected in the LDA. We will show in section HI that the 
LDA does give a spin stiffness, although formally its origin appears a bit differeB|p-| 

We now proceed to review the LDA approach of Kiibler, Sticht and co-workerg2al23. We begin by writing the ground 
state energy as a functional of the single-particle density matrix for a two-dimensional system in a constant external 
magnetic field 'B — Bz: 

+ \j j n{v)vH{v-v')n{v')(frcPr' + E^dPo..'{v)]. (2) 

Here, Ts[p(j ^'(r)] is the kinetic energy functional of non- interacting electrons, which in our case includes the external 
magnetic field B. The particle density is n(r) = TrpCT,o-'(r), w_f/(r) is the Hartree potential, and E^c [P(T,o-'(r)] is the 
exchange-correlation energy, which depends parametrically on the magnetic field (for ease of notation we omit this 
parametric dependence). At each point r we calculate an SU{2) matrix U(r) which locally diagonalizes pa-^a'i^)- 

^Ua,air)pa,a'ir)U*,^p{r) ^ SaJjUair), (3) 

with na{r) the eigenvalues of P(j,o-'(r). We write U(r) in the standard form 

/ e^'^Wcos^ e-^'^Wsin^\ 
U(r) = ... . ^2 ) _^ , s efr) - (4) 

where (^(r) and 9{r) are the local aziniuthal and polar angles of the magnetization density relative to a fixed coordinate 
system. The requirement that U(r) diagonalizes pa^ai{v) then yields 

tan^(r) = -lp4) (5) 

and 

2[(RepT,i(r))' + (ImpT,i(r))'l'^' 
tan6l(r) = — ^ — i , (6) 



[pt,t(i") -Pia(r)] 

with Re(2) and Im(z) denoting the real and imaginary parts of a complex number z. Equation (|^) gives 
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n^ir) = cos2 ^ + iptaWe^^^--) sin0(r) + ip|,^(r)e-'^W sm0(r) + p^^^ sin' ^ 
nx(r) = PT,T sin' ^ " ^PT.iWe^'^^"^^ sm(?(r) - ip|,^(r)e-^^W sin0(r) + p^a cos 



2 

2 0ir) 



(7) 



We now have a representation in which Po-.cr' (r) is locally diagonal, so that in the LDA we only need to know the 
exchange-correlation energy E-^c{^,S.) as a function of total density and polarization, or, equivalently, n-^ and n^. By 
using the SU (2) transformation, and by expressing the single-particle density matrix in terms of occupied KS orbitals 
ipi^a{r) we can then write the KS equations in the following form. 



T,V^,T(r) +wo(r)^^,T(I■) + Aw(r) [cos0(r)V'^,T(I■) + sin0(r)e''^('"V»a(r) 



T,V»a(i-) + «o(r)^.a(i-) + Av(r) 



)e-*'^W^,,T(r)-cos0(r)^,,4(r 



for each single-particle spin doublet ipi{r) = (V'i,t(r), Here, 

1 



Voir) = Woxt(r) + vu{r) + - [v^c^ir) + Wxca(r)] > 



and 



with 



Aw(r) 



(8) 



(9) 



(10) 



(11) 



In the KS equations there is a coupling between up- and down-spin components, so that in general the KS orbitals 
are now two-component spinors with both up- and down-spin components. Note that in an angular momentum 
representation, the equations (1q) couple z-components of orbital angular momentum Lz +mTi, with m an integer, and 
spin angular momentum Sz = to [L^, Sz — This coupling will also provide the niechanism for spin-charge 

textures in QHE systems in the same way as the Hartree-Fock equations by Fertig et al. doo. In the LDA, we write 



(12) 



where exc[i^, is the exchange-correlation energy per particle in an infinite, homogeneous system of filling factor v 
and polarization C = ('^T ^ '^l)/"-: ^nd i^(r) = 27rZ^n(r) is the density expressed as local filling factor. Then 



Wxc,t(i") = 

■fxcalr) = 



^ + -(i-o4 

^ + ^ fl? 



(13) 



with the derivatives evaluated at = v{v) and ^ = ^(r), so that 



d d 
vq{y) = Uoxt(r) + WH(r) + [ve^c{v,i)] - ^^e,,c{v 

Af(r) = ^exc(i',0- 



(14) 



Equations (H) wit|b-iWo(r) and Ai;(r) given by Eq. (|lj) are the KS equations which result from the approach by Kiibler 



and co-workerg23'E3. Here they are written in a form appropriate for the QHE. We now make the extension to an 
ensemble DFT by introducing occupation numbers fi^a for the two states a in each spin-doublet i, and by taking 

p,,,,(r) = V/,,„^l Jr)^t5^,^^_^(r), (15) 



where So- projects out the cr-component of the spinor ipi^ai^). IN ordinary DFT the occupancies fi^a are zero or 
one. In our ensemble DFT calcjalations, we obtain fractional occupation numbers using a method of running averages 
described in our earlier workMta. 
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III. FORMAL RESULTS AND NUMERICAL APPROXIMATIONS 



The usefulness of the LDA equations (ph, and (15) uhimately depends on the availabihty of good approx- 



imations for exc('^, C)i the exchange^correlation energy per particle of a homogeneous system of filling factor v and 
polarization ^. We have previouslyli3 described our first attempt at constructing an approximate energy surface for 
the QHE. In this section, we will describe in some detail our work to improve on that approximation. In particular, 
we have extended the exchange-correlation energy to the range v > \ and also better incorporated electron-hole 
symmetry. _ 

We start, as in our earlier workllj, by approximating t^^^y^^ as 

exc(^^, = 4c(^) + + ^''^^)- (16) 

Here, exc('^) is a smooth interpolation formula for the ground state energy of polarized QHE systems, e^(i') gives the 
cusps at the FQHE fractions, and /(^) is an interpolation formula obtained from considering only exchange in two 
dimensions, 

and bty^jyu) is the difference in energy between the fully polarized and the completely unpplapzed = 0) system 
at filling factor v. Some values for this latter quantity can be obtained from the literatureEZl"E2l, and the value for 
<5exc(i^ = 1) will be fixed later in this section using the spin stiffipess. We then use a spline fit to tie all these values 
of 5^t^(y) together|-t« a continuous function. We found earlierEB that the interpolation formula for et^iXy) given by 
Fano and OrtolaniEj together with our approximation for t^iy) gave excellent agreement between our DFT approach 
and numerical diagonalizationsc^ for small spin polarized systems. However, this interpolation formula was given 
only for v < \. To extend it to > 1, we have performed extensive numerical diagorLaJizations for toroidal spin- 
polarized systems of eight, nine, ten, and eleven particles in the two lowest Landau levels^. The data obtained from 
these calculations reveal a cusp in the exchange-correlation energy at v = 1. This cusp is due to the fact that the 
exchange-correlation energy per particle added to the second Landau level, after the lowest Landau level has filled up, 
is different from the exchange-correlation energy per particle for the filled lowest Landau level. We have confirmed 
this with analytical calculations (below), and constructed a simple analytical model to fit the numerical data. Figure 
^ depicts the data from the numerical diagonalizations and the analytical fit. 

We now briefly present the analytical calculation of the cusp in the exchange-correlation energy at i> = 1 for spin- 
polarized systems. We will do the calculation in a truncated Hilbert space, and take the state at = 1 to be a Slater 
determinant of the lowest Landau level single-particle wavcfunctions. This is the exact ground state of the system 
restricted to the lowest Landau level, and the interaction energy is Eq = Nciexdi' = 1,^ = 1) = ^^ci^/f^^, where 
A*"ei is the number of electrons. For iVoi -I- 1 particle we consider only eigenstates with N^i particles occupying a filled 
lowest Landau level plus one particle in the n = 1 Landau level. The lowest-energy state consisting of A'ci -I- 1 particles 
is then a linear combination of these degenerate Slater determinants with a uniform density. The direct (Hartree) 
energy of this state is canceled by a uniform positive background charge density. We can calculate the exact energy of 
the lowest-lying state (exact in the reduced Hilbert space used here) by considering the exchange interaction between 
a single particle in Landau level n = 1 at momentum fco = and a particle of momentum k in the lowest Landau 
level n — 0. This energy is (using the Landau gauge) 

ei(fc) = - / j d rd^r ——p^e ^'b e '"^e b ^—L — i\/(r - r') 

Hi{x/Ib) -^^^ ikv' -riri^'-^k? 
X — -e B eye b . (18) 

Here, Ly is the length of the system (taken to infinity at the end of the calculation), iJ„ is the n*'* Hermite polynomial, 
V{r) = e^ I (eo|r— r'l) the Coulomb interaction, and — l%k is the centerpoint of the harmonic oscillator of momentum 
k. The integrals in Eq. ( p^ ) can be evaluated to give the result 

eiik) = ^-^e-'^'^B/^^ [Ko{k%/4) + K,ikHl/4)] , (19) 

where Kn is the modified Bessel function of order n. We then finally integrate over all states k in the lowest Landau 
level and obtain 
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The cusp in ground state energy at v = 1 comes from the fact that the exchange energy of adding one particle to a 
system with a fuU lowest Landau level, ei, is different from the exchange energy per particle in the lowest Landau 
level. The cusp gives rise to a discontinuity in the chemical potential at = 1, which we need to evaluate. The 
chemical potential ai v — \~ is — yJiTple^ / [ttj, b) , and the chemical potential at = 1+ is 



ON,: 



= E{N,i + 1) - E{N,i) 



= [N,^eUv = 1, ^ = 1) + ei] - N,,e^,{v = 1,^ = 1) = ei = (21) 

2 V o eo«B 

Therefore, the discontinuity in chemical potential at = 1,^ = 1 is 

3 

Since the exchange-correlation potential for a polarized system is V^c(f^) — d [i^exc(i^, C — 1)] li'^^) = m(^) this discon- 
tinuity also appears in the exchange-correlation potential at v — \. 

In order to construct an exchange-correlation energy surface exc(i^, C) which gives a workable approximation for 
Q < V < 2 and < ^ < 1 we first construct an analytic approximation for the exchange-correlation energy at > 1 
and ^ = 1. We use a simple model in which we write 

exc(i^ > 1, e = 1) = a'^^e^^iu - 1, ^ = 1) + - + &. (23) 

V V 

This model is motivated by the fact that, as a first approximation, the interaction energy of a system with a full 
lowest Landau level and A^i particles in the n = 1 Landau level is approximately equal to the interaction energy of 
the full lowest Landau level plus the interaction energy of A'^i electrons in the lowest Landau level. The constants h 
and c then fixes the slope and value of the exchange-correlation energy at = 1+ and ^ = 1, and a is used to adjust 
this model to the numerical data. Fixing slope and value gives 

h=\e^,(y=\,i=\), (24) 

and 

c=-^exc(i^- 1,^=1), (25) 

respectively. A good fit to the numerical data is given by a = 2. 

Finally, we consider the exhange-correlation energy at ^ = 0. When the system is restricted to the lowest Landau 
level, particle-hole symmetry yields 

vt^,{v, i) - vt^,{v = 1, ^ = 1) = (2 - v)t^, [2 - V, {v^ - u^)/{2 - v)] - (2 - v)tUv = 1, ^ = 1), (26) 

where 

v = v^+vi>l (27) 

and 

e = ^^^^. (28) 

When the restriction of the system to the lowest Landau level is lifted, this symmetry is no longer exact because of 
inter-Landau level quantum fluctuations. However, we assume that it is only slightly violated, and construct Se^cW) 
so that exc(i^, respects this symmetry at ^ = 0. Using the form Eq. ( |l6|) and Eq. (^) at ^ = and v > 1 then 
yields 
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[2 - 1^] [el,{2 - + 6^(2 - ^) + Jexc(2 - i^)] - 2(1 - iy)e^,{,^ = 1,^ = 1) 

V 



(29) 



for V > \. Furthermore, at ^ = the exchange-correlation energy per particle has a continuous derivative at v = 1, 
which gives 



d((5exc(i^)) 



dv 



= -^yic{v = 1,^ = 1) - 5€^^{v = 1). 



(30) 



i/=i 



Equations (|T^), ([l7|), (p3|), (|29|), and (^0|), together with the data points for (5exc(i^) for i/ < 1 then define our 
exchange-correlation energy surface. A F90 subroutine package which evaluates the exchange-correlation energy and 
the exchange-correlation potentials for given (t/|, vi) is available from the authors. In this package, used in our 
calculations, we dropped the terms from Eq. (pa). This is of no consequence in our calculations presented here in 
which the total filling factor was never much greater than unity. 

To conclude this section, we show that the LDA does indeed give a spin stiffness; requiring this stiffness to be 
correct helps constrain de^di^ — !)■ We start by considering the total exchange-correlation energy of the system in 
the LDA, 



(31) 



where i^{r) and ^(r) are obtained from the local eigenvalues n-\{r) and ni{r) of the single-particle density matrixE3. We 
calculate in the LDA the change in exchange-correlation energy of an initially infinite, homogeneous, fully polarized 
system at = 1 in response to a gentle spin twist. The applied spin twist changes the local eigenvalues n^(r) of the 
single-particle density matrix, and thus iy(r) and ^(r). We write 



iy{r) = 1 + Siy{r), e(r) = l + '5e(r). 
The change in exhange-correlation energy is then (with n(r) = i^(r)/(27r/^)) 



1 

2^ 



Sn{r) 



[l + 6iy{r),l + S^{r)] (fr 



1 

2^ 



,(1,1) 



(32) 



(33) 



From the work of Moon et al. we know that the density variation is of second order in the gradient of the spin- 
rotation angle. Furthermore, general symmetry considerations give 



5i(r) cx [Vn{Y)Y , 



(34) 



where 0,{t) = z x m(r) is the angle through which the spin density is rotated, and m(r) is a unit vector parallel to 
the local spin density. We then expand A_Exc in powers of |Vr2|, and obtain to second order in |Vr2| 



AE^c = J <5n(r)exc(l,l)d'r + ^ J S,^{r) 



de 



dv 



9e„ 



(35) 



i/=i,j=i 



The point here is that an LDA in v{y) and ^(r) contains spatially varying polarization ^(r), and this corresponds to 
gradients in the magnetization density according to Eq. (|3J). The first two terms on the left-hand side describe a 
change in i?xc due a change in the density. For the spin twist we are considering 5n(r) (and 5v{v)) integrates to zero 
since no net charge is added to the system. The last term in Eq. (^5|) describes the change in E^^ due to a change in 
the polarization. This term gives the spin stiffness. In general, the spin stiffness ps is defined by the leading term in 
an expansion of the energy in gradients of the magnetization angle fl: 



Es 



1 



We obtain from Eqs. (M) and (p|) 



AE^r. oc — 



27r 



[wn{r)f d^ 



(36) 



(37) 
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By comparing Eqs. ( pq ) and ( |37| ) we see that ps oc {de^c/dS,), and it remains to work out the constant of proportion- 
ahty. To this end, we consider a system of spin-polarized electrons confined to the lowest Landau level. We use the 
Landau gauge A(r) = (0, Bx, 0) in which the single-particle basis functions are 



^(x — Xl.)'^ 



where 



Xk = l%k, k = 2TTn/Ly, n = 0,±l,±2, , 



In the initial state, all single-particle spinors are 



We take the spinors of the spin-rotated state to be 



i'k{x,y) 




cos%7/ife(a;,y) 



(38) 



(39) 



(40) 



(41) 



where A/c = 2'K/S.n/ Ly^ with An a fixed integer. This state is the rectangular analogue of a rotationally symmetric 
skyrmion. In order to calculate the exchange-correlation energy of this state within the LDA, we then need to find 
the eigenvalues of the single-particle density matrix 



Pa^a' (r) = 



Efc sin -f cos -fi^l^x, y)tpk-Ak {x, y) 
The eigenvalues are then readily obtained as 



X;fe cos2 ^|V'fe(a;, y)^ Y.k sin T ^-fi'kix, 2/)V'fe-A/c(a;, v) 

2 " 



Y^^siu' '-f\^k-Ak{x,y)\' 



(42) 



^cos^|0^,(x)+5:sin2|0. 



-Afc 



±- 



1 



^ cos^ y (/)fe(a;) - ^ sin^ -^(j^k-Ak 



2 ({fc, 

2 



k k 

To continue, we will make use of the following results 



sin 9k cos 6'fc(/ifc {x)4>k-Ak {x) 



1/2 



(43) 



2^11' 



and also the Taylor expansion 

(l)k{x)(f)k-Ak{x) ~ 4'k(.x) 



2ttI% 2 47r' 



0.(x)-Afc^^-H-(Afc) 



(44) 



(45) 



The second term on the right-hand side of Eq. (Bq) then vanishes when integrated over k. The third term contains 



(Afc)^=(|^) (An)^ 



which vanishes in the thermodynamic limit for fixed Art ^ 0. Therefore, we can approximate 

^(t)k{x)(l)k-Ak{x) K^(t)l{x) = - ^ 

k k 



2ttI%' 



(46) 



(47) 
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Next, we assume that 6k is slowly varying as a function of k. This means that in expressions like 



COS —(pu{x) — — ^ 
2 ^''^ ^ 27r 



dk- 



COS' 



''0{k) 



4>k{x) 



(48) 



where 4>k{x) as a function of k is sharply peaked about k = x/l^, we can expand to the trigonometric function as a 
function of k to second order in k about k — x/lg. When integrated over fc, all first-order terms containing dO/dk 
then vanish. Using Eq. (44) and d^6/dk^ — I'^d^ 6 / [dx"^) we obtain, after a little algebra, 



_ 1 1 
'1^'^ " 2 2^ 



1 

± - 



1 



2 27r?| 



^_il(dr 

4 \dx 



From the definition of the polarization ^(r) we then finally have 

2rii(r) 



<5^(r) 



1% fde 



4 V dx 



Inserting this into Eq. (BSh yields 



47r 



dx . 



(49) 



(50) 



(51) 



By comparing Eq. (36) and Eq. (|5l| ) we then obtain for the spin stiffness 



1 9ex 



47r 9^ 

In fact, the calculation is readily generalized to an arbitrary spin-polarized filling factor uq, so in general we have 



(52) 



Ps = 



(53) 



In our approximation for the exchange-correlation energy, we can then fix dexci^— 1) (for which there is no known 
value) by requiring that the LDA spin stiffness Eq. (p2h equals the known valueEj for the spin stiffness a,i v = 1, 



Ps 



1 1 

16 V27r eo^B 



This yields 



5exc{v = 1) 



7r2- ^/2 



0.1224 



(54) 



(55) 



As an indicator how good this value is for 5e^c{i^ — 1), we consider the state at = 1, ^ = in an approximation 
in which the spin-up and spin-down particles are completely uncorrelated. In this case the energy per particle of 
that state is the same as the ground state energy per particle at = i,^ = 1, which is aboutcj —0.469 e^/(eo'B)- 
Neglecting the correlations between spin-up and spin-down particles should lead to an overestimate of the energy per 
particle. This approximation yields a difference in energy per particle between v = = Q and = 1,^ = 1 of 
5tKc{v = 1) = 0.1577 e^/(eo^B)- In view of the fact that this is most likely an overestimate of the energy per particle 
at = 1, ^ = 0, we can conclude that the value 5e^c{v = 1) = 0.1224e^/(eo^B) obtained from the spin stiffness is very 
reasonable. The value 5e^c{v = 1) = 0.1224 e^/(eo/B) is the one we used in our calculations. 



IV. EDGE RECONSTRUCTIONS OF QHE DOTS 



We have applied the ensemble spin DFT approach to circularly symmetric QHE dots and studied the edge recon- 
structions of such systems as a function of edge stiffness. We use a model in which the confining potential is supplied 
by a uniform positive background charge density. The edge is modeled by a 'graded edge' in which the positive 
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background charge density goes to zero linearly as a function of radial coordinate r over a distance w. The total 
integrated positive charge is fixed and equal to the total electron charge; this then determines the radial distance 
at which the positive charge density starts to decrease. Even though concerns have been raised that this particular 
confinement is non-gcncrial3 it is the one which has been studied the most, and we chose it as a model confinement 
for comparisons with other work. _ 

It is known from Hartree-Fock calculationdl3 of = 1 Hall bars with a similar graded edge that for a sharp edge {w 
small enough), the electron gas is completely polarized and its density falls to zero abruptly near the edge. However, 
as the edge confinement softens, the electron gas develops an instability to .a spin textured edge for smaller values 
of g (the ratio of Zeeman to Coulomb energy), or a charge density waval2Hi3 for larger values of g, with the density 
modulated along the edge. Both instabilities break the translational invariance along the edge: in the spin-textured 
edge the spin density is modulated along the edge while the total charge density is constant; in the charge-density 
wave edge the spin density is constant along the_edge while the total charge density is modulated. For Hall bars at 
bulk filling v = 1/3, an effective field calculationE3 (here, Hartree-Fock calculations are obviously not applicable) also 
shows an instability to a spin textured edge as the confinement is softened. 

We would expect the analogous instabilities to occur for the circular dots. For stiff confinements, the electron 
density forms a so-called maximum density droplet (MDD), in which the electron gas is completely polarized with a 
filling factor which is unity in the bulk and which rapidly falls to zero at the edge (the MDD is the minimum angular 
momentum state for a spin-polarized system in the lowest Landau level). As the confinement is softened, the_edge 
should develop a spin textured or charge density wave instability. Indeed, our previous spin DFT calculationsEj, in 
which the magnetization direction was constrained to be fixed parallel to the external magnetic field, revealed that 
the MDD becomes unstable towards the formation of a partially polarized edge as the confinement was softened. 
This gave a variational bound (within the LDA) showing that the spin-polarized edge is not the ground state when 
the confinement is soft. The phase boundaries of the MDD that we obtained were in good qualitative agreement 
with Hartree-Fock calculations for Hall baralJ, although the obtained values of g at which the polarization was 
destroyed were much smaller than those from the Hartree-Fock calculation for Hall bars. We speculated that this 
difference is due to the different geometries or edge cqafinements. Our results were in rather good agreement with 
numerical diagonalizations using parabolic confinement^. The calculations presented here support the argument that 
the differences were due to the different confinements used. 

We have now extended our ensemble DFT calculations to include spin textured edges both for the MDD and for 
J/ = 1/3 droplets. Our new results show that when the edge becomes partially polarized, a spin textured edge has 
lower energy than one with constant direction of the spin density. In the calculations presented here, we have only 
considered states which do not break cylindrical symmetry of the charge density, and for which the azimuthal angle of 
the spin direction changes at most by 27r along any simple closed path. This excludes charge density wave instabilities 
and spin textures with topological charge greater than ±e, but imposing these symmetries simplifies the calculations 
a great deal. For example, the Hartree potential is easy to calculate for circularly symmetric charge densities, but 
considerably more tedious if that symmetry is broken. Nevertheless, the calculations we have performed satisfy the 
two most important criteria we wanted to establish: to demonstrate the usefulness of the spin ensemble DFT approach 
to general QHE fillings, and to establish through a variational bound (within the LDA) that for softer edges, the 
spin-polarized edge or an edge with constant spin quantization axis has higher energy than a spin-textured system in 
which the spin quantization axis is tumbling. 

We start with the ensemble spin KS equations Eq. (H) with the potentials vo{r) and Ati(r) given by Eqs. (|lj), and 
the single-particle density matrix found self-consistently using Eq. (^5|). We now make the simplifying assumption 
that the polar angle 9{r) of the spin density is a function of the radial coordinate alone, 0{r) = d{r), and that the 
azimuthal angle ipir) is of the form 

ip{r) = w0(r), (56) 

where v is an integer and 0(r) is the azimuthal angle of r in a planar polar coordinate system. In other words, we will 
restrict the modulation of the spin density along the edge to have a single Fourier component along the edge. With 



Eq. (56) inserted into the expressions Eq. (^ for vo{r) and At;(r) the spin-diagonal coupling conserves orbital angular 
momentum, while the spin off-diagonal coupling couples up-spin states ijjm+v,]- with angular momentum h(m + v) to 
down-spin states V'm,j, with angular momentum km: 



cos eirWmA (r) At;(r)e-™ (r) d'r. (57) 

This coupling is of the same form as that in the Hartree-Fock equations studied first by Fertig et oM. However, 
in the Hartree-Fock equations the off diagonal coupling was provided by an exchange integral, while our Eq. ( |57| ) 
also includes correlation effects within the LDA. Since the spin-diagonal coupling conserves angular momentum, the 
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diagonal elements of the single-particle density matrix are circularly symmetric so that the total density is circularly 
symmetric. We would like to point out that lifting the restriction Eq. ( p6|) allows for the possibility of breaking the 
circular symmetry of the charge density, which allows for charge-density waves in addition to spin-charge structures. 

We have solved the KS equations by expanding the spatial parts of each spinor V'i.Q(r) in the single-particle angular 
momentum basis functions 

^™,„(r) = L^-^L^^ fiL) e™^ (58) 

for the cylindrical gauge A(r) = ^B{xy — yx.) with LJ^ the associated Laguerre polynomials. We kept up to the n = 4 
Landau level in our calculations, and up to 120 angular momentum states (up to 1,200 single-particle states). For 
the V = 1 system we performed calculations of 40 and 70 particles in a magnetic field of 3.5 T with the bare Lande 
factor g* varying from 0.1 to 1.5 [g then varied from about 0.002 to 0.036, encompassing experimentally accessible 
values). The reason to keep the magnetic field relatively low, but still at an experimentally realistic value, was to 
fully include the effects of Landau level mixing typical in experimental system. We performed the calculations for 
V = 0,±1. The results can be summarized as follows: For small values of w, i.e. stiff confinement, the edge is spin 
polarized. As w increases, the v = Q channel becomes partially polarized. However, at the same value of w, the v ~ 1 
channel attains a lower energy with a nontrivial spin-charge texture. For the system sizes studied here, the value of 
w at which thej4nstability occurred is w ~ 11b- This is in quite good agreement with the Hartree-Fock calculations of 
Karlhede et aim. For a semi-infinite Hall bar, they found the onset to a charge-spin textured edge occurring at about 
w « 7Zb - 8Zb in the range of g from zero to about 0.03. As illustrations, figures || and ^ depict charge densities 
and spin rotation angle for w = 81b and g = 1.00. The v — channel (Fig. ^ has a small minority-spin density 
near the edge of the system, in qualitative agreement with our earlier calculations (the bump in density at the edge is 
characterisic of all confined QHE systems for which the confining potential is not macroscopically smooth). However, 
the V — 1 channel (Fig. which has lower energy, has a locally polarized spin density everywhere and a non-trivial 
spin texture, with the spin rotation angle 9{r) rising from to about tt/2 at the edge of the system. Note that for a 
bounded system there is no topological constraint on the spin rotation angle, as there is in an infinite system where 
the Pontryagin index has to be an integer. 

We have also studied a FQHE droplet at bulk filling of 1/3. While we are quite confident about our exchange- 
correlation energy for spin-polarized FQHE systems, it is not clear to us how good it is for arbitrary polarizations m 
the FQHE regime. For example, if we use the published value obtained from four-electron numerical diagonalizationsEZi 
for Se^c{i^ = 1/3) of 0.0017 e^/(eo's), the spin-stiffness in our model is J^155 x 10"'* e^/(eo'B), compared to the value 
of 9.23 X 10^^ /{cqIb) obtained from hypernetted-chain calculations^. It is not clear if the discrepancy between 
these two is mostly due to finite-size effects in the numerical diagonalizations, or our model Eqs. (|l|) and (|l^) of 
the exchange-correlation energy. We have also performed numerical diagonalizations of six electrons at v = 1/3. 
For six electrons, we obtain Se^cii^ = 1/3) = 0.00465 e^/(eo^B) - an increase of almost a factor of three from the 
four-electron result. This clearly shows that the energy per particle at total spin zero is much more sensitive to 
system size than the energy of the polarized ground state. Encouraged by our numerical results, we then we fixed 
<5exc('^ = 1/3) to 0.0136 e^/(eoZ_B) in order to have a simple model which gives a spin stiffness in agreement with the 
hypernetted chain calculations. We performed the calculation with 40 particles and the four lowest Landau levels, and 
up to 170 single-particle angular momentum states, both for the choice of confinement discussed earlier and parabolic 
confinement. 

Our calculations indicate that for a confinement provided by the positive background charge density, the system 
has an instability from a spin-polarized edge to a spiriptextured edge at an edge width oi w ~ 4:1b. Again, this 
compares rather well with the results of Karlhede et alx3 Using effective-field theories, they found an instability to 
a spin-textured edge at about w = 3.01b for g = 0.04. Note that, the effective-field theory tends to underestimate 
the value of w for which there is an onset to spin-textured edgest3. We also want to emphasize that in contrast to 
the effective-field (and Hartree-Fock) theory, our ensemble spin DFT is applicable to general inhomogeneous QHE 
system and includes the effects of Landau level mixings. As an example of our results, we show in figures ^ and |^ 
charge densities and spin rotation angles for a system with an edge width of w = 4 and Landc factor of g* = 1.00 in 
a magnetic field of 12 T (g « 0.04). The bump in total filling factor in the bulk of the v = channel (Fig. U) occurs 
quite generically for FQHE droplets, both for the choice of confinement discussed here and for parabolic confinement. 
In the case here, the system tries to take advantage of correlation energy at = 1/3 and = 2/5 to as large extent 
as possible as the edge is made wider. It does so by making the edge of the electron density sharper than that of the 
background charge, and making regions of = 1/3 larger than what is needed to accomodate all electron charge. The 
residual electron charge is piled up in a bump reaching v — 2/5. The v — 1 channel, on the other hand, is locally 
completely polarized with a spin rotation angle similar to the i> ~ 1 and w = & system for t; = 1. 
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We have also performed calculations for parabolically confined systems in the FQHE regime. The results are in 
this case more difficult to interpret with more complicated structures in the electron density and spin textures. For 
example, we have found that for some values of the magnetic field the v = channel can be a = 1/3 droplet with a 
large bump in density at the edge, while the v ^ 1 channel develops a hole with reversed spin density at the center 
and has lower energy for a small range of magnetic field. 

V. CONCLUSIONS AND SUMMARY 

We have here presented an ensemble spin DFT approach to general inhomogeneous QHE systems. This approach 
includes the spin degree of freedom, including non-coUinear spins, as well as Landau level mixing. On a formal 
level, we have demonstrated that our ensemble spin DFT in the LDA can give the correct spin stiffness at = 1. 
We have performed model calculations for circularly symmetric QHE dots in the ipteser and FQHE regime. These 
calculations show, in agreement with Hartree-Fock and cffcctivc-ficld calculationgl3il3, that the polarized system 
develops an instability as the confinement is softened and that the spin-textured edge attains a lower energy than 
the spin-polarized one. We have not included in our calculations the possibilities of charge-density waves, which may 
occurtllE3 instead of the spin textured edges at larger values of g. Preliminary calculations of parabolic dots show a 
surprisingly rich structure in spin- and charge-densities. This indicates that quantum dots in the FQHE regime is a 
rich subject yet to be fully explored. 

We have spent a great deal of effort on improving our approximation for the exchange-correlation energy exc(i', 0- 
At the present, we are confident that we have a very good approximation for spin polarized systems, and a good 
approximation for aribitrary polarizations and « 1. We are less confident about our exchange-correlation energy 
for arbitrary polarizations in the FQHE regime. Work need still to be done to refine the exchange-correlation energy 
for general FQHE systems. However, we are confident that the approach itself is robust and accurate provided good 
approximations for the exchange-correlation energy exist. 
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FIG. 1. Exchange-correlation energy vs. filling factor for a spin-polarized QHE system. The diamonds are results from 
numerical diagonalizations in the two lowest Landau levels, and the solid line for > 1 is our analytical fit Eq. (]2^). 

FIG. 2. Local filling factors i^]-{r), vi{r), and v{r) = h'^ {r) + i'i(r) in the ?; = channel for a 70 electron dot. The magnetic 
field strength is 3.5 T, g* = 1.0, and the edge width is w — S Ib- Here, the edge region is only partially polarized with a small 
minority-spin density. 

FIG. 3. The local total filling factor and spin rotation angle 6(r-) for the v = 1 channel of the same parameters as in Fig. 
^ This channel has a lower energy than the i; = channel, and is locally completely polarized with a non-trivial spin texture. 
That is, at any point r the spin density is polarized, but the spin direction is changing with position. 

FIG. 4. The local filling factor for the v = Q channel of a 40-electron dot in the FQHE regime. The external magnetic field 
strength is 12 T, g* — 1.0, and w — 4 Ib- In this case, the electron system is polarized. 

FIG. 5. The v = 1 channel for the same system as in Fig. ^ This channel is also completely polarized, but has a non-trivial 
spin-texture and lower energy than the v = channel. 
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